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ABSTRACT 


The linear boundary layer solution (Chang, 1973b), 
wren ineluvades the temporal acceleration, is utilized 
BOD eawGlonme parameterization to study tropical lower tro- 
pospheric waves. A one-level equation is derived for 
the top of the boundary layer and is solved numerically 
for its eigenvalues. Calculations are made for varying 
boundary-layer depths and vertical heating gradients. 
Growth rates, periods and eigenfunctions for the lowest 
meridional modes of synoptic-scale disturbances are 
analyzed. Results indicate that CISK heating does not 


prefer any finite wavelength. 
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i. ENT RODUCT ION 


Observational evidence of tropical lower tropospheric 
waves has been extensively collected and studied during 
the past several years. Although there is still some doubt 
moet Oo the detailed characteristics of these disturbances, 
it is generally observed that these waves propagate west- 
wea, exhibit a period of about four to five days in their 
meridional wind component, and may fall in one of the two 
wavelength ranges: around 2000-4000 km, or approximately 
8000 km. As a direct result of the observations, meteoro- 
logists have attempted to explain theoretically the 
mechanism and structure of these tropical disturbances. 

Three different hypotheses have thus far been postu- 
lated: 1) lateral forcing by mid-latitude systems, 2) 
meet kLOpac instability, and 3) thermal forcing. A concise 
but excellent description of these theories is found in 
feemmeada (1971). To summarize his diseussion briefly, various 
investigators have ee that all three hypotheses may be 
responsible for tropical wave growth. 

Widespread attention has been given the thermal forcing 
mechanism during the past decade. Ever since the inception 
Of the so-called CISK (Conditional Inston + ty Gf eine 
Second Kind) theory, first proposed by Charney and Eliassen 


(1964) and Ooyama (1964) to explain hurricane growth, many 





mavestigators have sought to extend its applicability to 
the larger scale tropical waves. Basically their CISK 
models incorporate thermal forcing by parameterizing the 
latent heat release from small scale convection in terms 
of the boundary layer convergence of the large scale flows. 
fire addition, Charney and Eliassen's geostrophic/Ekman 
beundary layer formulation (geostrophic vorticity and ver- 
mecal velocity are proportional at the top of the boundary 
layer) was commonly used. 

The primary approach in many of the previous studies 
Pero tO perform linear stability analysis of a system of 
equations in order to determine a preferred wave scale - 
that scale which has the largest growth rate when other 
Pemattions remain constant. Most of these earlier linear 
Studies used either a geostrophic/Ekman boundary layer or 
a Surface drag representation for the lowest level in the 
moecel. Although the linearization process was done to 
Simplify an analysis, a question generally arose in the 
Ekman-type boundary layer models near the equator. ino stlaas 
meagton the Coriolis pores in the equations of motion 
approaches zero and therefore the geostrophic relation may 
mo Longer be valid. re soe pointed out by Holton et al 
(1971) that a transition zone exists near the equator where 
mae Ekman boundary layer structure of ene more poleward 
Matitudes changes rather drastically into a Stokes-type 


boundary layer in the equatorial region. Generally, this 





mMeamsitton 2zOne Has not been accounted for by most of the 
investigators. 

For some of the earlier CISK wave models, a major dis- 
Grepancy exists. Yamasaki (1971) and Murakami (1972) 
found that preferred wave scales on the order of several 
thousand kilometers exist in their models. However, in 
@oemerast to these results, Chang (1971) and Williams and 
Robertson (1973) showed that a distinct maximum growth rate 
for the observed tropical waves did not occur. 

Because of this discrepancy, it remained inconclusive 
as to the effect CISK-type heating had upon tropical waves 
in the lower troposphere. The purpose of this research 1s 
to use an improved boundary layer representation within the 
peope Of linear theory to study the effect of CISK upon 
maeplcal wave growth. The linear boundary layer solution 
(Chang, 1973b) which includes the temporal acceleration is 
meea Lor the CISK parameterization. By neglecting the small 
temperature fluctuations in the thermodynamic equation, 
meen 1S generally valid for weak synoptic-scale disturbances 
in the lower tropical atmosphere, a one-level equation is 
derived for the top of the boundary layer. This equation is 
Similar to the free atmosphere equation in Chang (1973a). 
An eigenvalue problem is then formulated and the growth rates 


of the lowest meridional modes are obtained numerically. 





II. THE MODEL 


A. BASIC EQUATIONS 


The basic equations of the model are the linearized 


Bovssinesq primitive equations for an equatorial beta- 


plane: 











where 


ae ip 2 Ge _ 


u* Corry, 24 te 
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Zomale Mmerliadional;,- vertical and. time 
coordinates, respectively 


ZOnal” velocity 

Neri ronal weloci ty 
Vertical ve locivey 

pressure divided by density 


Meomuteal Eddy Viscosity scoefiicient, 
assumed constant 


local “change Of Coriolis parameter, F, 
WiHthiterespect to qdatertude, VY, at the 
equator; hence, equal to 22/a with & 
being the angular velocity of the earth 
and dae berng the radius -of the earth. 


In addition, the thermodynamic energy equation 


oT* 
ot* 
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* * 
mol + yee Se aes 9, (4) 


ox* dy* 





u 


10 





where 


T* = temperature 
ae _ g oT* 
Sewer tbactoe Stability parameter (= =— + ox ! 
, Z 
WiereowG ehGwtne acceleration Sis 


Mies GOsdmuavirty sand Gc is the 
Socelrleomwmeat Of dry alr at constant 
pressure) 
OF =-"dreabatic heating rate 
mo used. 
tine wlocal change and horizontal advection terms in (4) 
can be neglected considering the magnitude of each term 


mevaebace, 19/71; Holton, 1972). The thermodynamic equation 


then becomes 
wt*Sk* = Or (5) 


iis equation implies that the diabatic heating is balanced 
Meet he adiabatic vertical motion field. 

The basic equations are non-dimensionalized uSing a 

bp gr 
gravity wave speed, c = (gH) , where H is the scale height 


and the depth scale, and the Rossby parameter, $~. These two 


parameters are uséd to derive the non-dimensional units: 


- 1/2 
1) time, T = (cB) DEG Pimnorazontal length, L = €¢76) / ; 
| 2 
Speetiocrizontal velocity, c = L/t , 4) pressure, T7 = c /L , 
Setemperature, 9 = ae and 6) diabatic heating rate, 


3 
= Gi Rise 
W) Cy 


The non-dimensional equations are: 
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are assumed to be propagating waves 








Ap and tbequency™ vs 
Wey 2) 
V(y;,Z) 
ee et ASEVE) 405 
ply) 
Dae) 
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g oie ia) 
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2 = er 
Oz 
’ (3) 
(14) 


12 





Byeciltniaiatingeit,  ¥ and w in (11) = (14), a single 


Mercere lal —equcatrOn in p results: 


ap _ ,—#¥ _) 8p , A (y+ Dee 2) 
2 2 - 2° dy v 2 2 s 
dy y Vv (ya) 
2 2 
Ce ee le) 
a Vv OZ (3) : 2 


Mretnis equation the only vertically differentiated term 
is the heating term on the right-hand side. ras Hine lee 
that the system of equations is vertically coupled only 
mareough Glabatic heating. Hence, if the vertical profile 
iO 1S Specified, (15) can be applied as an ordinary dif- 
ferential equation with y as the only independent variable 


at any level above the boundary layer. 


B. THE CISK EQUATION 
The CISK process is parameterized so that the diabatic 
Piraerng LUunCtI1ON, §0lz), 1S related to the vertical velocity, 


W agit the top of the boundary layer (Z_): 


in T 


O4Z) «= 1 Cz) ° Wa ’ (16) 


where (Zz) is a Paoworstronality function which specifies 

Mme vertical profile of heating. Notice that (16) also 

implies negative heating when downward motion occurs at 

mame top Of the boundary layer. This is difficult to avoid 

due to the separation of the x-dependence of the solutions. 
The vertical motion at the top of the boundary layer 


is expressed by 


ies 





where 


and 
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and 


Coe od) 


2). i(v=y) 


Bear tiOon (17) is derived by integrating the total boundary 
Mmeeer COnvergence, uSing the solutions in Chang (1973b). 
This boundary layer solution includes the temporal accelera- 
mmom and, therefore, is frequency dependent. If the bound- 
ary layer depth becomes infinite (Z,, > O) ee eSNG U War dey ah 
PeeeewistsS at the critical latitude where y = v (Holton et al, 


T 


oy). iIf Zo fieatinate, the solution may possess quasi- 


Srmmgular behavior around this latitude depending upon Zn 
moomene mode of Symmetry about the equator (Chang, 19735). 
ine heating profile, n(z), is specified such that odn/doz. 
mimes from 1S to 2S at the top of the boundary layer. Con- 
Sequently (15) can then be applied at Zine [Notice that 
Mbeyeeand (16) require that n=S at the top of the boundary 
layer. For a lapse rate of 6°C eee and Wa cm oon the 
Magnitude of n in turn produces a heating rate of approxi- 
mately 3°C gay Because the slope of the heating func- 
tion is defined, it is implied that the heating function 
Maries smoothly in the vicinity of z = Ze Furthermore, a 


positive slope implies a heat source above the top of the 


bowndary layer. Finally, to reconcile the previous 


15 





assumption that the pressure perturbation is independent 


Sremerght for z°< z it is assumed that the smaller but 


mr! 
finite heating in the uppermost layer of the boundary 
layer acts 1n such a way that p remains constant in height. 


A second-order homogeneous differential equation in p 


is cbtained by Substituting CLG) eandevOl/ wel nto tl oa : 


E ee OE oe + G p = QO 


: (18) 


where E, F, and G are functions of y and \Vv: 


2 2 
2 Cy dn 
= Se —t 
4 VS al az 
2 2y ye) R a 
2 2 eS) 2 age 
(7 9 a 
A (yee UN) 2 i) an 
B24 y eae Y = Ee 
{y - Vv ) 
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Pe. So Go. ON PROCEDURE 


tomo om Ls jemumectcaliyv, a central difference scheme 
in terms of the independent variable, y, is utilized. The 


difference expression is 


ee) itis, (ay) 7-2 By py (BF, 54) py = 0 29) 
Ue peo y Seu ea 
where 0 LS the total number of grid points and 
Rep width between the two boundaries 


y=! 


feenmannel—-£Llow boundary condition is assumed such that the 


meeturbatlons vanish at the boundaries; that is, 
eee (20) 


The Fee eae equation, (19), is solved for its eigen- 
Mabue, V, using an iterative procedure as described by 
Stone TAG TORS Biesote bore a given) value of 7A; Z mnt K and 
dn/fdz, an initial guess for v is made. Then by satisfying 
Wieowouethe DOundary conditions in (20), the pressure per- 
migeatronuat the first Grid point is set to zero; that 1s, 
Py = 0. Next, the amplitude of the eigenfunction at the 


Zecomemoriad MOLnt 1S arbitrarily chosen as unity. Using 


17 





these first two values for the eigenfunctions, the numeri- 
cal solution process successively solves for the remaining 
values of the eigenfunctions: Por Ee 
In general, the guess for V is incorrect and the calcu- 
lated value for De is therefore non-zero, thereby violating 
mame second boundary condition of (20). Consequently, the 
eateculation of Ps 1S repeated, using a value of V which dif- 
fers from the initial guess by 1%. The two values of Pp, are 
then used to perform a linear extrapolation from the two 
Peamues Of VY to a new value of V which would make I) 220) If 
it depends on Vv ina linear fashion. The new value of Vv is 
then used to reiterate the whole procedure, and iterations 


are repeated until a value of V accurate to 0.01% is ob- 


tained. In summary the process is as follows: 


(2) 


BVA oe P 


ge Sales eas . 


The upper limit for N is prespecified as 25 in almost all 
mye computer runs. Convergence to a solution generally 
mBecurs well before this upper bound but, if convergence is 
not achieved by N=25, that particular run is terminated 


and a new set of input values is automatically provided to 


Maxdimize Computer efficiency. 
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Two external parameters are varied in the following 
range: A= 1 km to 4 km and odn/dz = 18 to 2S. The solu- 
tions are obtained for both the symmetric (pressure is 


symmetric about the equator) and the asymmetric (pressure 


is asymmetric about the equator) cases. 


1B, 








Ve RESULTS 


A. FREE MODE 

Solutions for the "free" mode, when the right-hand 
emecre, Of (iS) as set to zero, are obtained as a test for 
ihe model. Zero growth rates (infinite e~folding times) 
result for all wavelengths. 

In Figure 1, the results for wave period are shown. 
Moetice that wavelengths of the synoptic scale range (2000 
~10,000 km) have periods between approximately three and 
seven days for the asymmetric case and five and ten days 
for the symmetric case. In general, the longer waves 
have higher frequencies than the shorter waves in both 
cases. These free mode frequencies correspond quite 
erosely to those described by the baretropic Rossby wave 


frequency formula 


eee 


where 2 = 1/D, D being the width of the channel, and the 


approximated Matsuno (1966) wave-frequency formula 


v= ry (Ro 2n e+ «lL ) & 


Here k is the non-dimensional zonal wave number and n is 


the meridional mode number. 
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'. GROWTH RATE 

Preures 2 through > indicate the growth rates for both 
Symmetry cases. THe wimeportant result shown by all four of 
these figures is that the synoptic scale waves have essen- 
tially the same growth rate. Thus, there exists no pre- 
ferred wavelength in the synoptic scale range - a result 
which agrees with those of Chang (1971) and Robertson (1973). 

For each set of external parameters, two different solu- 
tions are generally found: a long-wave (LW) mode, one that 
Pessesses greater growth rates in the long-wave portion of 
Ehe Synoptic scale; and a short-wave (SW) mode, one that 
favors short-wave growth. Both modes of the symmetric case 
generally show larger growth rates (faster e-folding times) 
than those of the asymmetric case, although for shorter 
wavelengths this difference rapidly decreases. 

Preeiaerease in the slope of the heating function, 1; 
increases the growth rates of the waves. For a boundary 
layer depth of one kilometer, both modes of the asymmetric 
and symmetric cases indicate an approximate doubling of 
their growth rates with an increase of dn/dz from 1s to 2g 
(compare emnee 2x Ana: <4.) Box Zon = 2 km, the same increase 
Mmieche vertical gradient of heating produces an even larger 
change in the growth rates - generally a 125% increase 
results (see Figures 3 and 5). 

Likewise, a deeper boundary layer produces larger growth 


rates. This effect is especially noticeable for the LW 


ran) 





mode of the asymmetric case. Grewth cates for this. mode 
approximately triple when the boundary layer deepens from 
one to two kilometers, but increase more slowly with 
deeper boundary layers (results not shown). The SW mode 
for the asymmetric case, as well as the two symmetric 

case modes, exhibits similar but less pronounced increases 
in growth rates when the boundary-layer depth increases 
from one to two kilometers. 

A rather interesting result is that the symmetric SW 
mode for shallow boundary-layer depths are Soi Oana 
wavelength of 4500 km), a low vertical heating gradient 
(dn/do2 = 1S) and longer (>4000 km) wavelengths has larger 
growth rates than the LW mode (Figure 4). But for deeper 
imoundary layers the opposite occurs; that is, the growth 
Bate Of the symmetric LW mode is greater than those of the 
SW mode (Figure 5). Thus, a simple relationship between 
the growth rates of the symmetric modes and the boundary- 
iavyer depth 1s not apparent. On the other hand, Figures 
2 and 3 indicate the asymmetric SW mode has larger growth 
rates than those of the LW mode for both boundary-layer 


depths. 


c. PERIOD 

AS was mentivened earlier the periods of the free modes 
(Figure 1) were determined as a test case. Figures 6 
Porough 9 allustrate the effects of heating and boundary- 


layer depth upon wave period. Note that the periods for 


Ze 
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the LW modes of both the asymmetric case and symmetric 
mocciomuvery schosely resembple those of the corresponding 

free modes. This result is expected, because the LW modes 
generally have much smaller growth rates than the SW modes, 
thereby better approximating the free modes' zero growth 
rates. Both the LW and SW modes of the asymmetric case 
exhibit shorter periods than those of the corresponding 
modes of the symmetric case. 

It is also noted that varying on/dze produces little 
change in frequency. Only a slight decrease in period is 
noticeable with an increase in odn/odz. The deeper boundary 
layer does not appreciably affect the wave period either, 
although a very small increase in frequency is detectable 
for the longer wavelengths of the symmetric LW mode and 
a small decrease for the asymmetric LW and SW modes, as 


well as for the symmetric SW mode. 


1 ELGENFUNCTIONS 

Bagures 0 .and 11, which represent the free mode pro- 
moles for ae =i and 2 km, respectively, are consistent 
Pen the solutions in Chang (1973b). When CISK heating 
is included, the p and Win Meotiles Lome both symmetry cases 
generally tilt northeast-southwest between the northern 
boundary and the equator. However, the equatorward limit 
morethe tilt Hn the peprofile of the symmetric case depends 


upon wavelength (the longer the wavelength the farther away 


aN) 





meom the equator is the limit). Such a northeast-southwest 
mec in the eigenfunctions implies poleward transport of 
westerly momentum and agrees with observations summarized 
meewallace (1971). 

Somparison Oot wFigures 2 and 16 (among others) shows 
that the tilt for the asymmetric LW mode decreases with 
longer waves. The SW mode generally shows a similar change. 

By increasing the vertical heating gradient at the top 
of the boundary layer, the tilt increases for both asymme- 
tric modes especially in the short waves (see Figures 12- 
ik). The only significant change in tilt due to an increase 
in the boundary-layer depth occurs for the LW mode, for 
which the tilt increases with decreasing wavelength (compare 
Pees J2 with 14 and 13 with 15). 

Figures 12 and 16 show a decrease in the SW mode's 
egquatorward maximum of Win amplitude and a general smoothing 
of the we profile when the wavelength decreases. On the 
other hand, the LW mode shows a slight equatorward shift 
of the ue Maximum and a uniform amplitude increase with 
meereasing wavelength. The p profiles for both modes re- 
Main nearly unchanged in these figures. 

An increase in odn/dz shifts the asymmetric LW mode's 
maximum in ores and p more equatorward, particularly for the 
Shorter waves (See Figures 12 and 13). The SW mode, mean- 
while, shows a slight equatorward shift only in its pole- 


ward We Maximum and a somewhat smoother profile of ire when 
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dn/dgZ increases. The p profile of the SW mode remains al- 
mest unchanged with a larger dn/dz. 

Deepening the boundary layer also effects an equator- 
ward shift in the asymmetric LW mode's p and W on maxima 
Geempare Figures 12 with 14 and 16 with 18). ahs ee sul & 
fmemesimilar to that found by Chang (1973b). From Figures 
16 and 18, it 1S apparent that for longer waves the SW 
mode's equatorward maximum in Wi becomes the primary maxi- 
mum when the boundary-layer depth increases. While both 
ae Maxima of the SW mode shift slightly toward the equator 
with a deeper boundary layer, the p profile remains 
feeetudlly unchanged (Figures 16 and 18). Both modes indi- 
cate a larger ao amplitude when om increases from 1 km to 
ee while the p amplitude again remains nearly constant 
(eempare Figures 12 with 14 and 16 with 18). 

Essentially the same results as those just noted for 
the asymmetric case are found for the symmetric case. 

Only the significant differences are described below. 

From Figures 24.and 25 it is seemmunat the egquatorward 
ee maximum of the symmetric SW mode shifts poleward for 
the longer wavelengths when on/oz increases. in comtras te, 
the asymmetric SW mode's equatorward a maximum remains 
stationary when dn/dozZ is increased but the poleward maxi- 
mum moves equatorward. While an increase in the boundary- 
layer depth shifts the asymmetric LW mode's W on and p maxima 


more toward the equator, only the w = maximum moves 
mn 


At 
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equatorward for the symmetric LW mode (compare Figure 20 
Ween 22 and 24 with 26). 

The phase difference between the pressure and vertical 
motion fields generally remains 1/2 cycle for both symme- 
try cases under the conditions investigated. Tire Jon ly 
exception appears with the symmetric case. Pimehnes VY lel way 
of the critical latitude and equatorward, a rapid phase 
Sart in Wp, occurs such that the two fields are nearly in 
Phase at the equator. Whereas an out-of-phase relation- 
ship implies a downward energy transfer from the source 
region immediately above the boundary layer, in-phase im- 
pimees the converse transport of energy, that is, from the 
boundary layer to the atmospheric layer immediately above. 
However, the relative amplitudes of both p and Me are 


small near the equator, thereby Minimizing the upward 


transport of energy within the near-equatorial region. 
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Ve SUMMARY SAND CONCLUSTONS 


A linear, tropical boundary-layer model is formulated 
On an equatorial beta-plane by utilizing the primitive 
equations of motion, a simplified thermodynamic energy 
Equation and a CISK parameterization. Solutions in the 
monrm Of complex wave frequencies are obtained numerically 
aS eigenvalues to a second-order finite difference equa- 
tion. 

The primary conclusion reached is that there is no 
Preferred synoptic-scale wavelength for CISK-forced tropi- 
cal waves. Although two modes for each symmetry case 
exist when the CISK forcing iS incorporated. in no case 
does a noticeable peak in the growth rate curves 
Materialize in the synoptic scale range. tn aed tien, 
meme e@-folding times for the more realistic vertical heat- 
mig gradient (dn/oz = 18) are found to be of the order of 
5 days or slower, even for the deeper boundary-layer 
depths. This magnitude of growth rate is quite small 
compared to baroclinically unstable waves in the midlati- 
tudes. Consequently, even if one scale is slightly pre- 
ferred among the others, its relatively faster growth 
would be at such a slow rate that the GompetbTvoLon, Lor 
energy among the waves may not be strong enough to pro- 


duce a dominant wavelength. 
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iiewictinmeote lUStOn Of this study - that the growth 
rate is nearly uniform in the synoptic scale range - agrees 
with some earlier works (Chang, 1971,and Williams and 
Robertson, 1973), but differs from others (Yamasaki, 1969, 
ime7is> Hayashi, 197l;and Murakami, 1972). A similar inde- 
mendence of scale 1s noted for the so-called "spin-down" 
Pmeoolem. by Specitying dNi/dz = -1 and substituting it into 
mee Lorecing f£Unction on themrightaehand side of (15), a 
crude formulation of the spin-down problem is made. By so 
@oung, the vertical motion profile for the free atmosphere 
1S assumed to decrease linearly from le at the top of the 
boundary layer to zero at one scale height (1.e., z2 = H). 
The resultant damping rates are found by using the same 
peocedure for computing the CISK growth rates and they are 
shown in Figure 28. The spin-down times are small as 
expected because of the small rotation rates near the 
equator. ineaacaltion the curves in Figure #ze. show nearly 
uniform spin-down rates for the synoptic-scale waves. The 
CISK model equation, therefore, is mathematically analogous 
to the spin-down equation except that the forcing function 
has the opposite sign. Hence, it may not be surprising 
Ghat the CISK growth rates are also independent of scale. 

Although there is no preferred wave scale in the 2000- 
10,000 km range, the LW mode of both symmetry cases never- 
theless exhibits larger growth rates in the synoptic scale. 
PMimeentnmasct, the short wave end of the LW mode has decreas- 


ing growth rates with decreasing wavelength. This mode 
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with its growth rate characteristics seems to agree with 
the HB mode of Yamasaki (1969) and the findings of Murakami 
(1972). Pmt tion che ir results are Similar to the LW 
mode in the following two ways: 1) maximum amplitude 
eeeeotc tear the equator, and 2) wave periods of about 5 
ways OCCUr. Onwenieceotner hand, “the SW mode, especially 
for the shallower boundary layers, has its maximum ampli- 
Eude more poleward and possesses larger periods. The 
pmemode also appears to prefer zero scale for growth. 
These three characteristics indicate that the SW mode may 
memraentified as a hurricane (or typhoon) mode. 

From the results of this study, a growth rate depen- 
dence upon boundary-layer depth appears to be two-fold. 
Mig@csctee by increasing the depth of the boundary layer, 
Seowem rates increase for both the LW and SW modes of both 
symmetry cases. Secondly, the greatest change in the 
growth rates for both cases generally occurs between Z 
= 1 km and oe = 2 km. As a result of these findings, it 
Mmemapparent that linear, tropical boundary-layer models 
must account for the variations in boundary-layer depth 
mrieeorder Eo, clenemae more accurately the growth rates of 
Em@Oopical lower tropospheric waves under CISK forcing. 

Within the tropical boundary layer, there is generally 
a transition zone between a quasi-~Ekman boundary layer of 
the more poleward latitudes and a guasi-Stokes boundary 


layer of latitudes near the equator. The slow-rotation 
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Stokes boundary layer with its inefficient pumping does not 
@enmeribute much to wave growth. Ene winclisron of the tem= 
poral acceleration terms in the boundary-layer equations 
motos Seems tO have very little effect in CISK instability. 
fmeeract the results of this study generally agree with the 
gGuasi-~geostrophic models of Chang (1971) and Williams and 
Robertson (1973). 

At the top of the boundary layer, the vertical motion, 
Be assumes a smoother profile across the transition zone 
Wemitical latitude) for the growing solution compared to 
miagemrree SOlution in Chang (1973b). However, the large 
Pmese Shift across the transition zone remains Similar to 
mmae £Lound for the free solution. 

iwimernedar adveckive terms may be important in “atfect— 
ifamgeehe Growth rates of Synoptic-scale tropical waves. 
Such terms are neglected in this study and should be in- 
emided in future research to determine their effect ina 
Gjeeok model. Finally, a possible improvement to this linear 
model may be the inclusion of lateral shear of the mean 


Meow tm the equations for both the boundary layer and the 


free atmosphere. 
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Figure 1. Wave periods for free modes. 
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Solid lines indicate pressure perturbation, 
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Eigenfunctions for 4500 km wavelength, 
asymmetric case, = 2 km and on/dz = 2s. 
Solid Jines eat pressure perturbation, 
Py Telemcashed lines the vertical motion, 
Won at the pop. or Cie boundary layer. 
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Eigenfunctions for 2000 km wavelength, 
symmetric case, 2,= 1 km and on/oz = 1S. 
Solid lines indicate pressure perturbation, 
bp, ana cached lines the vertical motion, 
Went at the top of the boundary iayer. 
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Eigenfunctions for 2000 km wavelength, 
symmetric case, Z2_= 1 km and on/oz = 2S. 
Solid lines indicate pressure perturbation, 
p, and dashed lines the Verticalemetion, 
Went at the top of the boundary layer. 
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Eigenfunctions for 2000 km wavelength, 
symmetric case, zZ,= 2 km and On Zan eos 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical Metron, 

Mire at the top of the boundary layer. 
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Eigenfunctions for 2000 km wavelength, 
symmetric case, z,= 2 km and On 22 ao 
Solid lines indicate pressure pertrubation, 
p, and dashed lines the vertical motion, 
ea at the top of the boundary layer. 
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Eigenfunctions for 4500 km wavelength, 
symmetric case, 2,= 1 km and ane, aoe Ge 
Sold lines andicate pressure perturbation, 
p, and dashed lines the vertical motion, 
Wo Abe tne toprer tile boundary Jayexr. 


54 





LATITUDE 


megure 25. 





AMPLITUDE PHASE 


Bigentunc elens fOr 4500) kml wavelength, 
symmetric case, z,= 1 km and ony doze = 22. 
Sotid lines: indicate pressure perturbation, 
p, and dashed lines the Vertlcal Mot lon, 
Wit at the top of the boundary layer. 
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Eigenfunctions for 4500 km wavelength, 
symmetric case, Z,= 2 km and dn/dz = 1s. 
Solid lines indicate pressure perturbation, 
Dy aud dashed slimes the vertical motion, 


w_, at the top of the boundary layer. 
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Eigenfunctions for 4500 km wavelength, 
symmetric case, 2,= 2 km and gn/foz = 2S. 
Solid lines indicate pressure perturbation, 
pb, and dashed lines the vertical motion, 
w_, at the top of the boundary layer. 
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